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Abstract. We compare the history of the galaxy mass build-up, as inferred from 
near-IR observations, and the Star Formation Rate of massive stars in the comov- 
ing volume traced by deep extensive far-IR surveys, both possible now with the 
Spitzer Space Telescope. These two independent and complementary approaches 
to the history of galaxy formation consistently indicate that a wide interval of cos- 
mic epochs between z ~ 0.7 to ~ 2 brackets the main evolutionary phases. The 
rate of the integrated galaxy mass growth indicated by the IR-based comoving 
SFR appears consistent with the observed decrease of the stellar mass densities 
with redshift. There are also indications that the evolution with z of the total pop- 
ulation depends on galaxy mass, being stronger for moderate-mass, but almost 
absent up to z = 1.4 for high-mass galaxies, thus confirming previous evidence 
for a "downsizing" effect in galaxy formation. The most massive galaxies appear 
already mostly in place by z ~ 1. 

Although a precise matching of this galaxy build-up with the growth of nuclear 
super-massive black-holes is not possible with the present data (due to difficulties 
for an accurate census of the obscured AGN phenomenon), some preliminary 
indications reveal a similar mass/luminosity dependence for AGN evolution as for 
the hosting galaxies. 
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1. Introduction 

The cosmological origin of the galaxy mor- 
phological sequence can now be very effec- 
tively constrained by the wide IR multi- 
wavelength coverage offered by the Spitzer 
Space Telescope, including the near-IR 
photometric imaging by the IRAC and the 
far-IR imaging by the MIPS instruments. 
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By combining these with the unique imag- 
ing capabilities of HST / ACS and the enor- 
mous photon-collecting power of spectro- 
graphs on large ground-based telescopes 
(VLT, Keck), we have now a definite chance 
to directly picture in-situ the process of 
galaxy formation. 

The latter has remained a rather con- 
troversial issue untill recently. Published 
results from high redshift galaxy surveys 
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appear not unfrcqucntly in disagreement 
with each other (see Faber et al. 2005 
for a recent review). This is partly due 
to the small sampled areas and the corre- 
sponding substantial field-to-field variance. 
However, a more general problem stems 
from the apparent conflict between reports 
of the detection of massive galaxies at very 
high redshifts (e.g. Cimatti et al. 2004; 
Glazebrook et al. 2004) and indications for 
a fast decline in the comoving number den- 
sity at z > 1 (e.g. Franceschini et al. 1998; 
Fontana et al. 2004). 

The most direct way of constraining 
the evolutionary history of galaxies and 
trying to resolve these discrepancies is to 
derive the redshift-dependent stellar mass 
functions from deep unbiased surveys (e.g. 
Bundy, Ellis & Conselice 2005). We con- 
tribute to this effort by exploiting in this 
paper very deep public imaging by the 
GOODS project to select a complete sam- 
ple of high-z (z < 2) galaxies selected in the 
3.6 /urn IRAC near-IR band. Near-IR sur- 
veys are best suited for the study of faint 
high-redshift galaxy populations, for vari- 
ous reasons. Compared to UV-optical se- 
lection, the observed fluxes are minimally 
affected by dust extinction. At the same 
time they are good indicators of the stellar 
mass content of galaxies (Dickinson et al. 
2003; Berta et al. 2004), and closer to pro- 
vide a mass-selection tool. For typical spec- 
tra of evolved galaxies, the IRAC channel 
also benefits by a K-correction particularly 
favourable for the detection of high-redshift 
galaxies. 

In addition to the Spitzer observations, 
the GOODS and related CDFS projects 
have provided the community with an un- 
precedented amount of high quality optical 
and near-IR data in CDFS, particularly the 
very deep multi-band ACS imaging, allow- 
ing the most accurate morphological anal- 
ysis currently possible. We will illustrate 
here the power of combining such multi- 
wavelength information in the analysis of 
the evolutionary mass and luminosity func- 
tions of faint galaxies. 



A complementary view on galaxy for- 
mation is possible by direct sampling the 
rate of star formation based on suitable 
tracers. As established by exploratory ob- 
servations with ISO and SCUBA (e.g. 
Franceschini et al. 2001; Elbaz et al. 2002; 
Smail et al. 2002), most accurate SFR de- 
terminations require the detection of dust 
emission in the far-IR, which often includes 
the majority of the radiant energy by mas- 
sive stars. Our direct estimates of the evo- 
lutionary stellar mass functions arc then 
compared with the most updated results 
on the comoving density of star forma- 
tion by Spitzer far-IR observations. Modulo 
an assumption about the stellar initial 
mass function, the two approaches provide 
the complementary integral and differential 
views of star formation, respectively. 

Galaxy formation has involved a vari- 
ety of complex phenomena (see e.g. Baugh 
et al. 2005), including tidal interactions 
and merging, black-hole formation, accre- 
tion, and feedback processes. The ubiqui- 
tous presence of super-massive black-holes 
in galaxy nuclei (Magorrian et al. 1998) 
and the relevant observed scaling laws 
(Ferrarese and Merritt, 2000) suggest that 
feedback by nuclear activity (Springel et al. 
2005) might have had an important influ- 
ence on the process. This motivates us in 
looking also for comparison with the evo- 
lutionary patterns in the Active Galactic 
Nuclei population. 

We adopt a standard cosmol- 
ogy (r^M=0.3, f^A = 0.7) and ex- 
press the H dependence in terms of 
h = H o /W0 Km/s/Mpc. 



2. Cosmic Evolution of the Galaxy 
Mass and Luminosity Functions 

Of the two approaches to the history of 
galaxy formation - the one based on the 
study of the instantaneous SFR with suit- 
able indicators, and the other estimating 
the evolutionary stellar mass functions - 
the latter has been recently recognized to 
benefit by various advantages (e.g. much 
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Fig. 1. Image at 3.6 fim of the GOODS region in the CDFS. The area covered is 12 x 18 square 
arcminutes approximately and the exposure time per sky pointing was 23 hours as a minimum. 
The image sensitivity is better than 1 \ijy for point sources. 5622 sources with Sz.n^ > lfJ-Jy 
are detected inside the Spitzer/ACS common area of 160 sq. arcmin, 5302 of which are galaxies 
and 320 stars. 



lower extinction uncertainties) and to be 
more robust. 

The evolution of the galaxy stellar 
mass function has been recently studied 
by us using a multi-wavelength dataset in 
the Chandra Deep Field South (CDFS) 
area, obtained from GOODS (Dickinson et 
al. 2004) and other projects (VVDS, Le 
Fevre et al. 2004) , and including very deep 
high-resolution imaging by HST/ACS. Our 
reference catalogue of faint high-redshift 
galaxies, which we have thoroughly tested 
for completeness and reliability, comes from 
a deep (53.6 > 1 A*Jy) image by IRAC 
on the Spitzer Observatory (see Figure 
|y. These data are complemented with ex- 
tensive optical spectroscopy by the ESO 
FORS2 and VIMOS spectrographs, while 
deep K-band VLT/ISAAC imaging is also 
used to derive further complementary sta- 
tistical constraints and to assist the source 
identification and SED analysis. We have 
selected a highly reliable IRAC 3.6/^m 
sub-sample of 1478 galaxies with S3. 6 > 



10 //Jy, 47% of which have spectroscopic 
redshift, while for the remaining objects 
both COMBO- 17 and Hyper z are used to 
estimate the photometric redshift. 

Based on this extensive dataset, we 
have estimated time-dependent galaxy lu- 
minosity and stellar mass functions, while 
luminosity/density evolution is further con- 
strained with the number counts and red- 
shift distributions. The deep ACS imaging 
allows us to differentiate these evolutionary 
paths by morphological type, which can be 
reliably performed at least up to z ~ 1.5 
for the two main early- (E/S0) and late- 
type (Sp/Irr) classes. 

These data, as well as our direct es- 
timate of the stellar mass function above 
Ahh 2 = 10 10 M Q for the spheroidal sub- 
class, consistently evidence a progressive 
dearth of such objects to occur starting at 
z ~ 0.7, paralleled by an increase in lumi- 
nosity. A similar trend, with a more mod- 
est decrease of the mass function, is also 
shared by spiral galaxies, while the irregu- 
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Fig. 2. Mass (right-hand panels) and luminosity (left-hand panels) function estimates derived 
from the 3.6 /J,m IRAC/GOODS sample with S3. 6 >10 fj,3y, splitted into three redshift bins from 
z — 0.1 up to z = 1.4. The contributions of the various morphological classes is marked with 
different symbols: early-types (open diamonds - dotted lines), late-types (open squares - dashed 
lines), total (filled circles - solid lines). The thin solid line on the right marks the local mass 
function from Cole et al. (2001). 



lars/mergers show a positive evolution (in- 
creased number) up to z ~ 1.5. 

It is interesting to note that the de- 
crease of the comoving density with red- 
shift of the total population appears to de- 
pend on galaxy mass, being stronger for 
moderate-mass, but almost absent until 
z = 1.4 for high-mass galaxies, thus con- 
firming previous evidence for a "downsiz- 
ing" effect in galaxy formation (e.g. Cowie 
et al. 1996). 



based on optical observations of the UV 
rest-frame flux (Lilly et al. 1995; Madau et 
al. 1996; Le Fevre et al. 2004). However, it 
has become recently evident that a prob- 
ably major fraction of the emission by 
the most massive, luminous and short-lived 
stars, during the stage when they are still 
embedded inside their parent dusty molec- 
ular clouds, is optically extinguished and 
reprocessed at wavelengths from the mid- 
IR to the sub- millimeter. 



3. Probing the Evolution of the 
Comoving Star Formation Rate 
with Far-IR Observations 

The usual way of determining the star for- 
mation activity in high-rcdshift galaxies is 



3.1. The ISO and SCUBA explorations 

The Infrared Space Observatory mission 
and the SCUBA/JCMT observations have 
offered the first tools for systematic cos- 
mological surveys at long wavelengths. 
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22 'x22 ' field 

total time: 43 ksec 

redundance factor ~7 

283 sources detected 
between 0. 1 and 8. 
mjy (Ssigma) 

9% of sources are 
stars 

85% of optical 
identifications in r' 
(r'<24.5) 



Fig. 3. Image at 15 /im of an area of 22x22 arcmins in the Lockman Hole by Rodighiero et al. 
(2004). Details about the observations are reported in the figure. 



The good imaging capabilites of ISO have 
been exploited for deep mid-IR observa- 
tions with sensitivities sufficient to detect 
sources at cosmological redshifts (Elbaz et 
al. 1999, 2002; Franceschini et al. 2001). 
Surveys in the wide band at 12-18/im 
(\ e ff = 15 /zm) have been performed in 
the ISOCAM GT, over a total area of 
1.5 square degrees, with >1000 sources de- 
tected (Elbaz et al. 1999). Rodighiero et 
al. (2004) and Fadda et al. (2004) report a 
deep and a shallow survey of the Lockman 
Hole. An image of the deep field is reported 
in FigureOl The two Hubble Deep Field ar- 
eas (North and South, total of ~ 50 sq. ar- 
cmin) have been deeply surveyed at 15 fim 
down to 100 fiJy (Rowan-Robinson et al. 
1997; Aussel et al. 1998). At brighter fluxes, 
ELAIS has observed a total of 12 square de- 
grees at 15 \im (Oliver et al. 2000; Vaccari 
et al. 2005). 

The differential 15/im counts of extra- 
galactic sources revealed a sudden upturn 
at S15 < 3 mjy, and a later convergence 
below Si 5 - 0.3 mjy (Elbaz et al. 1999; 
Gruppioni et al. 2003). These counts, far in 
excess of the no-evolution predictions, re- 



A 2-phase SF: origing of galactic disks and spheroids 

Quiescent star-formation Active bursts of SF 

• Quiescent phase of SF • Transient starburst phases, 
during most of the recurrently triggered by 
Hubble time, galaxy mergers and 

• slowly building stars interactions 
from the regular flow of • Violent relaxation, 
gas in rotational generation of galaxy 
supported galaxy disks spheroids 

• Thin-disk geometry => • Large-scale redistribution of 
moderate dust the dusty ISM into the inner 
extinction, originates the galactic regions => high 
optical BKG extinction => origin of CIRB 

Fig. 4. A two-phase scheme for the evolution 
of galaxies selected in the optical and the IR. 



quired strong evolution of the galaxy lu- 
minosity functions. At the faintest mid-IR 
fluxes, the redshift distributions of the ISO 
galaxy population showed an excess num- 
ber of luminous sources between z=0.5 and 
z=l. The upper redshift boundary is due to 
the ISOCAM mid-IR camera responsivity, 
implying a strong K-correction penalty for 
z>l. 
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Fig. 5. Panel a: comparison of the differential 15 |im galaxy counts measured by ISO and 
Spitzer (Teplitz et al. 2005), showing excellent agreement. Panel b: relative contribution of low 
luminosity starbursts (Lib, < 10 11 Lq), IR luminous sources (LIRGs, Lib. > 10 11 ) and ULIRGs 
(Lir < 10 11 ) to the total SFR density of the universe as a function of redshift (Perez-Gonzales 
et al. 2005). Panel c: evolution of the comoving IR galaxy emissivity up to z=l (green region) 
and the respective contributions from low luminosity galaxies (Lib < 10 11 , blue), IR luminous 
sources (Lib > 10 , orange) and ULIRGs (Lib < 10 11 , red). 



These results have been analyzed by 
Franceschini et al. (2001) and Elbaz et al. 
(2002) in conjunction with data and lim- 
its on the integrated sky intensity between 
1 and 1000 /zm set by COBE observations 
(the CIRB background) , and with the sub- 
mm galaxy counts at 850 /im by SCUBA 
(see review by Smail et al. 2002). The com- 
bined multi-wavelength constraints imply 
that the far-IR volume emissivity of galax- 
ies increases strongly from local to red- 
shift z ~ 1 and should flatten or converge 
above. A comparison of the results of opti- 
cal and IR cosmological surveys has been 
interpreted by Franceschini et al. (2001) 
to indicate that galaxies during their life 
form stars both in a quiescent mode, mostly 



responsible for the optical emission, and 
through short-lived starbursting episodes 
originating the far-IR emission, as summa- 
rized in Figure 0] 



3.2. Results by the Spitzer Space 
Telescope 

The far-IR imagers (MIPS and IRS) on 
Spitzer not only have confirmed the ISO 
results (Teplitz et al. 2005, see panel (a) in 
Figure EJ, but, thanks to the longer wave- 
length capabilities of the 24 /xm channel, 
have been able to extend the ISO surveys 
to higher redshifts cover the critical redshift 
interval of z ~ 1 to ~ 2. 
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Fig. 6. Evolution of the comoving IR luminos- 
ity density and SFR as a function of redshift, 
taken from Perez-Gonzales et al. (2005). All 
values in the figure are computed for h — 0.7. 




Fig. 7. The red-shaded area shows the ex- 
pected evolution of the total integrated stel- 
lar mass density, as inferred from measure- 
ments of the comoving SFR. The upper dat- 
apoints and the continuous line are direct esti- 
mates by Franceschini et al. (2006) for the total 
galaxy population. The lower (cyan) points are 
the corresponding measurement for spheroidal 
galaxies, showing a somewhat faster conver- 
gence with z. 



A region of the CDFS was observed 
at 24 by Le Floc'h et al. (2005) dur- 
ing the MIPS instrument GTO over a to- 
tal area of 0.6 deg 2 and PSF FWHM of 6 
arcsec. A sample of 2600 sources brighter 
than 80 pJy was combined with existing 
optical data in the field and used to de- 
rive bolometric IR luminosity functions and 
SFR's from z=0 to ~ 1. These results im- 
ply that the comoving IR energy density of 
the Universe evolves like (1 + ^) 3 9±0 - 4 up 
to z ~ 1 (panel (c) in Fig. |5j). 

From MIPS 24 pm observations of the 
CDFS and HDFN, complemented with a 
systematic photometric redshift analysis 
exploiting the Spitzer IRAC data, Perez- 
Gonzales et al. (2005) confirm on one side 
the fast increase of the comoving SFR den- 
sity to z = 0.8. As shown in Figure |£3 
the SFR density is found to continue ris- 
ing from z ~ 0.8 to 1.2 with a smaller 
slope, but then to remain roughly constant 
above, in quite good agreement with results 
by Franceschini et al. (2001). 

4. Discussion 

It is of interest to verify how consistent 
the results about the redshift evolution of 
the integrated comoving stellar mass den- 
sity pm(z) in Fig.|2]and the integrated star 
formation rate psfr(z) in Fig. [f)|are. Both 
quantities are based on the assumption of 
a universal Salpeter-like stellar IMF. The 
relationship of the two is simply: 

Pm{z) = Pm(0) - / psFB,(z)—dz (1) 
Jo dz 

where psfr.(z) ~ 0.014 h (1 + z) 3 9 at 
z < 1 and constant above (from Perez- 
Gonzales et al.), and where we can approx- 
imate dt/dz = l/[H (l + z) 2 - 5 ]. Assuming 
p M (0) ^ (6.8 ± 0.8)10 8 h M Q /Mpc 3 from 
Fukugita et al. (1998), we derive the evo- 
lution of the integrated stellar density re- 
ported as red shaded region in Figure 
This is compared with the direct estimates 
by Franceschini et al. (2006) for the to- 
tal galaxy population (suitably scaled from 
their fig. 19 to correct for the different 
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mass boundary). Within the large uncer- 
tainties, we see here that the two comple- 
mentary views of the history of star forma- 
tion appear to provide consistent results, 
hence in broad agreement with the assumed 
Salpeter IMF. 

Panels (b) and (c) in Fig. 03 also re- 
veal some interesting features in the z- 
evolution of the SFR density as a function 
of the galaxy luminosity. Low-luminosity 
sources display their main SF activity be- 
low z = 1, while the highest-L ULIRGs de- 
velop their main activity phases at z > 1.5. 
Intermediate luminosity objects peak at in- 
termediate redshifts of 0.7 to 1.5. This, to- 
gether with the results illustrated in Fig.El 
provide concordant evidence that the for- 
mation and evolution of galaxies has been 
a strong function of luminosity and mass, 
as anticipated and discussed by various au- 
thors (Cowie et al. 1996, Franceschini et al. 
1998, 1999; Treu et al. 2005). 

Various attempts have been proposed 
to explain this effect, which seems to in- 
dicate a sort of inversion in the process of 
dynamical assembly of galaxies compared 
to the hierarchical expectation. One possi- 
bility considers energy feedback by a nu- 
clear AGN stopping gas accretion and star 
formation in the most massive systems at 
high-z (Granato et al. 2004; Springel et al. 

2005) , while leaving it undisturbed in the 
lower mass galaxies where AGN activity is 
irrelevant. 

This interpretation would call for a 
close relationship of the evolutionary his- 
tories of star-formation and nuclear Black 
Hole Accretion in AGNs. Unfortunately, a 
precise matching of galaxy build-up and 
the growth of nuclear super-massive black- 
holes is not possible with the present data, 
due to difficulties for an accurate census of 
the obscured AGN phenomenon. The prob- 
lem has not yet been finally settled by the 
inclusion of data by the Spitzer surveys, ba- 
cause even a detailed IR multi-wavelength 
coverage is not sufficient to identify the 
whole of the obscured AGN population 
(see Franceschini et al. 2005; Polletta et al. 

2006) . 



AGNs and SBs in SCUBA Galaxies 




10" 10" 10" 10" 10" 10" 
Rest-frame 0.5-8 keV luminosity (erg 3'') 



Alexander et al. 2005 



Fig. 8. Relation of far-IR to X-ray luminosity 
in higfi-z SCUBA galaxies (from Alexander et 
al. 2005). The bulk of these (red circles) are 
found to require both starburst and AGN si- 
multaneous emissions. 

Preliminary results, mostly from X-ray 
surveys, confirm on one side that AGN 
and starburst activity occur concomitant 
in luminous high-z forming galaxies (as 
shown e.g. in the remarkable outcome by 
Alexander et al. 2005 reported in Figure 

EJ. 

It is also interesting to note that a sim- 
ilar mass/luminosity dependence for AGN 
evolution as for the hosting galaxies (with 
fast evolution with cosmic time for the 
highest-L and deferred activity for the 
lower-L objects) has been found in the lu- 
minosity functions of type-1 AGNs derived 
from complete X-ray surveys by Hasinger 
et al. (2005). 
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